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Abstract 

Consider a random three-coordinate lattice of spherical topology having 2v vertices 
and being densely covered by a single closed, self-avoiding walk, i.e. being equipped 
with a Hamiltonian cycle. We determine the number of such objects as a function of v. 
Furthermore we express the partition function of the corresponding statistical model 
as an elliptic integral. 
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1 Introduction 



Hamiltonian cycles play an important role in the study of the configurational statistics 
of polymers P] . Given a lattice, a Hamiltonian cycle is defined as a closed, self-avoiding 
walk which visits each vertex of the lattice once and only once. Typically, in polymer 
physics, one has been interested in calculating the number of Hamiltonian cycles for a 
lattice with a certain regular structure and a given number of vertices. This counting 
problem has been exactly solved only in a few cases, namely for the two-dimensional 
Manhattan oriented square lattice 0, the two-dimensional ice lattice 0] and the two- 
dimensional hexagonal lattice 0, 0. A field theoretical approach to the problem of 
counting Hamiltonian cycles on regular lattices was invented by H. Orland et. al. 
and has recently been further developed by S. Higuchi |7|]. 

Here we address the problem of counting Hamiltonian cycles on random lattices. 
Whereas on a regular lattice the choice of boundary conditions has to be carefully 
considered and plays a crucial role for the outcome of the analysis (see for instance 
for a discussion) on a random lattice no such complications occur. Furthermore one 
might speculate if the complex dynamics of polymers does not call for a random lattice 
rather than a regular one. 

Consider a random three-coordinate lattice of spherical topology having 2v vertices 
and being equipped with a Hamiltonian cycle. In the present paper we aim at deter- 
mining the number of such objects as a function of v. The objects in question can 
also be viewed as spherical triangulations consisting of 2v triangles and being densely 
covered by a single closed, self-avoiding walk (where it is understood that the walk 
consists of links belonging to the dual lattice). Not any three-coordinate graph has 
a Hamiltonian cycle. Obviously, one-particle- reducible graphs do not, but there also 
exist one-particle-irreducible graphs which can not be equipped with a Hamiltonian 
cycle. We show an example of such a graph along with the corresponding triangulation 
in figure 1. In order to study the above mentioned combinatorial problem we consider a 
matrix model which generates graphs or triangulations of the relevant type, namely the 
n — ► limit of the fully packed 0(n) model on a random lattice. The 0(n) model is a 
so-called loop-gas model and the loops, being closed, self-avoiding and non-intersecting 
play the role of the walks. In the fully packed model all vertices of the generated graphs 
are visited by a loop and in the limit n — > only one loop remains. The regular lattice 
version of the fully packed 0{n) model was studied in references [||, ||, [5|, 

The 0(n) model on a random lattice comes in two different versions, an unrestricted 
one [|1(], [Tl[ and a restricted one . In the latter the loops are restricted to having even 
length while in the former there is no restriction on the loop length. On a regular lattice 
all loops necessarily have even length. The two different versions of the random lattice 
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(a) (b) (c) 

Figure 1: A three-coordinate one-particle-irreducible graph, which does not have a 
Hamiltonian cycle (a), and the corresponding triangulation (b) or equivalently (c). 



0(n) model have different critical properties. The restricted 0(n) model belongs to 
the universality class of the unrestricted model fll2| . The fully packed, restricted 
0(n) model with n = 2 is of particular interest because it describes the generalisation 
of Baxter's three-colour problem [JT3] to a random lattice |TJ|, |T2[]. As long as one 
considers triangulations consisting only of an even number of triangles (such as those 
corresponding to a two-dimensional manifold without boundaries) the n — > limit of 
the fully packed, unrestricted 0(n) model coincides with the n — > limit of the fully 
packed, restricted 0{n) model. Since the unrestricted 0{n) model on a random lattice 
is mathematically easier to handle than the restricted one, we shall take as our starting 
point the unrestricted one. We note, however, that our n — > solution applies to the 
restricted model as well. 



2 The Model 



Our starting point is the following matrix model 

Z(g)= [ dA 17 dBi exp I -N Tr 
Jnxn f = \ { 



-i -i n n 

i=l i=l 



(2.1) 



where A and Bi are hermitian N x N matrices. This matrix model generates tri- 
angulations which are densely covered by closed, self-avoiding and non-intersecting 
loops which come in n different colours. More precisely we can write the free energy 
F{g) = log Z(g) as 



F(g) = T,N- 2h F h (g), 

h=0 
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F h {g) = 52M h (2v)g\ N h {2v) = Y j n l M^{2v). (2.2) 

v=l i=0 

The quantity (2v ) counts the number of genus h triangulations built from 2v tri- 
angles and densely covered by i closed, self- avoiding and non-intersecting loops 0. In 
the following we shall aim at determining N'q 1 \2v), i.e. the number of spherical tri- 
angulations consisting of 2v triangles and being densely covered by a single closed, 
self- avoiding walk. First, let us write instead of (|2.1| ) 



2(9) 



(Vtf 2 z(g) 



[ dAT\ dBi exp I —N Tr 

Jnxn fJi { 



i -i n n 

—A 2 + - T]Bf — Y] AB 2 
2g 2 ^ 1 ^ 



i=l 



i=l 



(2.3) 



Furthermore, let us introduce a resolvant W(z) corresponding to the model given by 
the partition function Z (g) 



1 1 00 
W{z) = (- Tr -> = £ N~ 2h W h (z). 

We shall be interested only in the large- N limit of the resolvant, Wq(z) 



W (z] 



Z—>CG 



m=0 



yJTl+l 



Let us consider the quantity t 2 

2g 2 d 



U 



(^TrA m ), =0 . 



d 



(2.4) 



(2.5) 



log{(Vg) N2 Z(g)} h=o =g + 2g 2 -F (g) 



N 2 dg 

OO V 

9 + 2EE«% (,) (^)^^. 



(2.6) 



v=0 i=0 



We see that we can achieve our aim of calculating jVq (2t>) by determining the contri- 
bution to £2 which depends linearly on n. In order to do so, we shall make use of the 



well-known loop equation for the resolvant ||10| , pT| . This equation can be derived by 
making use of the invariance of the matrix integral (|2.3|) under analytic redefinitions 
of the fields and reads in the large- iV (spherical) limit 

- (zW(z) + (1 - z)W(l -z)-2) = (W(z)f + (W(l - z)f + nW(z)W(l - z) (2.7) 

g 

We now introduce the following functions 

1 



+ 



3)' 



+ 



z m (l - z)p (l-z) m zP' 



{21 



4 Here it is understood that a triangulation, T, with a loop configuration {£} is counted with the 
weight l/Aut(T{£}) where Aut(X{£).) is the order of the auto-morphism group of the triangulation T 
with the loop configuration {£} 
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The functions ^ m , p {z) can be expressed in terms of the £k(z) as follows 

UM = t ( m + P p Z \- k ) 6W + 1 ( ™ + l 1 1 " * ) & W (2-9) 

and in terms of the ^-functions we can write the loop equation fl2.7| ) as 

1 00 n 00 

- ^6 — ^2tkUik+i+2 + - X! ik+i,i+itkU- (2-io) 

# z=i jU ^ fcj=0 



Using the relation fl2.9|) and identifying the coefficients of each £ m in ( P-lOj ) we get for 

m > 2 

1 m ~ 2 00 f k + l \ 

-t m = £ tktm-k-2 + ft J! tm+h-ltl ( 7 ] • (2-H) 
# fc=0 Z,fc=0 V / 

This equation can be solved perturbatively in n and g. First, let us set t m = g m ^ 2 T m . 
Then we have 

m— 2 00 v 

fc - 1 



— £ TkT m -k-2 + ft £ £ ( l _ 1 ) T v -kT m +k-2 g v ^ 2 (2.12) 



k=0 v=l k=l 

From this we see that the TVs have the following expansion 

OO V 

T* = ££ft'TfV /2 (2-13) 

and that 

2vAf^(2v) = T^ 2v) . (2.14) 

3 Exact solution as n —> 

For n = our model (|2.3| ) is nothing but the Gaussian one-matrix model and we have 

T (o,o) ( 2k V- T (o,o) _ f3 ^ 

Here is the Catalan number of order A; which counts in particular the number of 
distinct configurations of k non-intersecting arches drawn on top of a line. 
Equating the terms linear in n in equation ( |2.12|) we see that 



k=l 



J2(2v-l)]_v ( v + 2 \( v -l\ 
(v + 2)\(v-l)\t 1 {k+l ) [k-l) 

2(2w - l)!(2w + 1)! 



(v + 2)\(v + l)\v\ (f-1)!' 
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To prove the last equality one can for instance compare the terms of power zero in x 
on the two sides of the obvious identity 

- (i + *r 2 - (i + -Y" = -^r(i + x?'" +l • (3.3) 

x x \ xj x v+L 

Next, let us define a partition function for spherical triangulations densely covered by 
a single closed and self-avoiding walk 



v=l 
oo 



v - (2r-l)!(2r+l)! v 

— V Wnfr+M \ g *». (3.4) 

1g 2 h (v\f(v + l)\ (v + 2)(v + l)v y y } 



Here g plays the role of a cosmological constant. Using |TI| we find 
9" fn^Wf^ _ l^(2«)!(2t; + l)!^ 



dg i v / 5 ^ (t>!f (t> + 1)! 

= J-_ [if (4^) - E{4y/g)] - - g for < g < 1 (3.5) 

where if and i? are the complete elliptic integrals of the first and second kind re- 
spectively. The radius of convergence of the series above is g = g c = A. Using the 
analyticity properties of the elliptic integrals one gets 

FP{g)~^— 2 (g c - g f\og{g c -g) for g ~ g c . (3.6) 

This shows that the value of the critical index 7 defined by F^^g) ~ — g) 2 ~ 7 takes 
the value —1 which is the value characteristic of c = — 2 conformal matter coupled to 
quantum gravity. This scaling is in accordance with the scaling obtained in |]L6[|. 



4 Direct combinatorial approach 

The number A/q (2i>) can be derived more directly from the following purely combina- 
torial argument. For each spherical three-coordinate lattice covered by a Hamiltonian 
cycle, let us choose one particular link visited by the cycle. For lattices with 2v ver- 
tices, there are exactly 2v such visited links. By cutting the link and pulling apart 
the two created end-points, we can deform continuously on the sphere the (now open) 
Hamiltonian cycle into a straight line. Since the path is Hamiltonian, all the vertices 
of the original lattice lie on this line. To fully re-construct the lattice configuration, we 
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simply need to specify the image under the deformation of all the unvisited links. Since 
we are dealing with a three-coordinate lattice, each vertex on the line is the extremity 
of exactly one such unvisited link. After the deformation, unvisited links thus form a 
system of v non-intersecting arches which connect the 2v vertices in pairs, see figure 2. 
Pairs can be connected either above or below the straight line. Note that by choosing 
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Figure 2: Left: A triangulation consisting of 18 triangles, equipped with a Hamiltonian 
cycle. Right: The corresponding systems of arches. 

which side of the line is the top or the bottom side, we artificially double the number 
of configurations. 

Under the above deformation, the same arch configuration is obtained exactly 
Aut(T{£}) times, where Aut(T{£}) is the order of the auto-morphism group of the 
original triangulation T with the loop configuration {£}. This factor cancels out the 
pre-factor l/Aut(T{£}) in the definition of Mq 1 \2v). All distinct arch configurations 
are thus counted exactly once. We can finally identify 2x2vx A/" (1) (2 v) as the number 
of distinct systems of v non-intersecting arches connecting 2v vertices on a straight 
oriented line from above and/or below this line. 

By decomposing the v arches into a system of k arches above the line and v — k 
arches below, we get 

4v^ 1 \2v) = jr[ 2 2 V k y k c v „ k , (4.1) 

where the binomial coefficient comes from the specification of the vertices which are 
connected from, say, above and where c k denotes the number of distinct configurations 
of the k arches above the line. This number is nothing but the Catalan number given 
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in (|3.1|), which leads to 



AvM^(2 




{2v)\{2v + 2)\ 
(v + 2)\((v + l)\) 2 v\ 



(4.2) 



in agreement with (|3.2|). 



5 Conclusion 

We have solved the non-trivial combinatorial problem of determining the number of 
spherical triangulations consisting of 2v triangles and being densely covered by a single 
self-avoiding and closed walk. Let us define the entropy exponent of such objects ujh, 
as 

logu H = Km i- log (Mi 1] (2v)) . (5.1) 



From ( |3.5| ) we see that 

lu h = 4. (5.2) 
For triangulations without decorations the corresponding exponent takes the value 

JP = 2 • 3 3 / 4 (5.3) 

and for triangulations corresponding to one-particle- irreducible three-coordinate graphs 

= VT (5 ' 4) 

As mentioned in the introduction, one-particle-reducible graphs do not contribute 
to Afi 1) (2v) since they do not have any Hamiltonian cycles. Not any one-particle- 
irreducible graph contributes to A/q (2v) either, but those which do are counted with 
a weight which equals the number of ways they can be equipped by a Hamiltonian 
cycle. 

We call the attention of the reader to the fact that we have obtained an exact 
expression for the spherical contribution to the free energy of the restricted as well as 
the unrestricted, fully packed 0(n) model on a random lattice in the limit n — > 0. An 
exact solution of the complete, restricted 0(n) model is still lacking. Such a solution 
would amongst others provide us with an exact solution of the three-colour problem on 
a random lattice [0, . The results of the present paper might be viewed as a step 
towards the exact solution of this problem. To finish, we may ask how important is the 
Hamiltonian nature of the walks, or more generally whether or not the fully packing of 
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Figure 3: Construction of the SOS degree of freedom from the fully packed loops, on a 
regular triangular lattice and on a random vertex three-colourable triangulation. The 
thick lines indicate links whose dual link is occupied by a loop. The height variable 
h, defined on the vertices of the lattice, is such that the difference of height is ±1 for 
nearest neighbours connected by a thick line and ±2 for nearest neighbours connected 
by a dashed line. The sign is fixed by demanding that h and the number i of the 
corresponding sub-lattice T, are equal modulo 3. For the regular lattice, this SOS height 
variable corresponds precisely to the underlying picture of a 2d piling of elementary 
cubes. 



the (restricted or unrestricted) 0(n) model is a relevant constraint. For an 0(n) model 
on the regular honeycomb lattice, it was explained in || that the fully packed fixed 
point differs from the low temperature dense phase fixed point by the emergence of an 
additional SOS degree of freedom which can be constructed if and only if the loops are 
fully packed. This degree of freedom is defined as a height variable h v living on the 
vertices v of the dual triangular lattice (see figure 3). Indeed, this lattice is naturally 
divided into three alternating sub-lattices 7$, i = 1,2,3. Fully packed loops can then 
be translated into height variables by demanding that h v = i mod 3 for v G Tj and by 
imposing for nearest neighbours v and v' that \h v — h v >\ — 1 if the link between v and v' 
has its dual link occupied by a loop and \h v — h v >\ = 2 if not. This extra SOS degree of 
freedom results in a shift by 1 of the central charge of the fully packed fixed point with 
respect to the usual fixed point for the dense phase. On a random lattice, vacancies 
(i.e. unvisited vertices) can be introduced for instance by adding a term proportional 
to A 3 in the potential in ( p.l|) . As shown in |Tl|, on a random lattice the nature of the 
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dense phase fixed point is not changed by eliminating vacancies and the constraint of 
fully packing is thus irrelevant in this sense. This is also what we find here for n — > 
where we get c = —2 for Hamiltonian walks on a random lattice, as for the dense phase 
of the 0(n — > 0) model and not c = — 1 as for the fully packed 0(n — > 0) fixed point 
describing Hamiltonian walks on the honeycomb lattice. This result is not surprising 
since, on a random lattice, fully packed loops cannot be transformed anymore into an 
SOS height variable and we thus expect no shift in the central charge. As a natural 
question, we may ask whether it is possible to preserve the SOS variable on a random 
triangular lattice. A SOS variable can be defined in the same way as above if the 
random lattice can be divided into three alternating sub-lattices. In other words, one 
needs that the vertices of the lattice can be coloured in three colours such that nearest 
neighbours have different colours or, stated in yet another way, that any given vertex is 
shared by an even number of triangles. A vertex three-colouring problem of the above 



mentioned type was discussed and solved in |18| . Its coupling to a fully packed 0(n) 
model is still an open question but, according to the above discussion we expect the 
same shift in the central charge as for the regular model. Note that fully packed loops 
automatically have even length on a vertex three-colourable triangular lattice. The 
case n = 2 is of particular interest since it would describe a problem of three-colouring 
of both the vertices and the edges of random triangulations, a problem equivalent to 
the folding of random triangulations, as explained in |12|, [18] . Based on the discussion 
above, we expect c = 2 for this model, as for the fully packed 0(2) model on a regular 
lattice §. 
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